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OBJECT  OF  TASK 

To  Improve  existing  knowledge  on  gamma  and  neutron  shielding 
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ABSTRACT 

An  analytical  approach  is  developed  to  permit  determination  of 
ganma  radiation  attenuation  as  it  passes  through  two-legged  rectangular 
ducts  and  shelter  entraneeways.  The  approach  used  employs  the  albedo 
concept  for  wall  scattering  and  includes  correction  terms  necessary  to 
account  for  the  "comer  lip  effect."  With  appropriate  simplifying 
assumptions,  moderately  simple  engineering  formulas  are  obtained. 

Actual  use  of  the  formulas  requires  better  knowledge  of  differential 
angular  albedo  than  Is  presently  available;  however,  by  assuming  Iso¬ 
tropic  distribution  of  the  albedo  function,  a  very  good  comparison  of 
experimental  Information  with  results  calculated  by  this  technique  Is 
obtained. 
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I.  INTRODUCTION  TO  THE  PROBLEM 

The  problem  of  designing  a  shield  against  nuclear  radiation  such  as 
gamma  rays  and  neutrons  la  essentially  one  of  Interposing  sufficient 
mass  of  material  between  the  source  of  radiation  and  the  detector. 

Both  gamma  rays  and  neutrons  exhibit  an  approximate  exponential  atten¬ 
uation  behavior  In  passing  through  matter,  so  that  complete  removal  of 
all  the  radiation  Is  impossible.  The  amount  of  radiation  has  to  be 
reduced  to  a  level  consistent  with  the  purpose  of  the  shield.  The 
design  of  such  a  shield  would  be  relatively  simple  If  either  the  source 
of  radiation  or  the  shield  space  could  be  completely  surrounded  by  the 
proper  thickness  of  shielding  material.  This  is  not  possible  since 
practical  conalderatlona  require  that  some  means  of  access  to  the 
source  or  shielded  space  must  be  made.  Thus  a  reactor  core  must  have 
control  rods,  conduits,  etc.  A  personnel  shelter  must  have  utility  and 
ventilation  ducts  and  personnel  entranceways. 

The  problem  la  to  be  able  to  design  these  access  ducts  so  that  the 
level  of  radiation  outside  Che  radiation  source  enclosure  or  inside  a 
shelter  Is  hot  significantly  Increased  by  the  installation  of  these 
ducta  and  openings.  Zt  has  frequently  In  the  past  been  Che  practice 
merely  to  place  ona  or  two  90-degra(i  bends  In  ducts  without  any  very 
careful  analysis,  to  provide  necessary  attenuation.  With  rather  thick 
shielding  walls  and  small  ducts  this  practice  has  proved  adequate. 
However,  for  a  large  duct,  such  as  an  entranceway  Into  a  shelter,  this 
crude  approach  is  Inadequate.  The  design  of  such  an  entranceway  will 
probably  depend  greatly  on  Its  nuclear  radiation  transmission  proper¬ 
ties,  and  since  entranceways «are  an  expensive  component  of  a  shelter, 
careful  analysis  Is  very  worthwhile. 

It  Is  Che  purpose  of  this  note  to  analyze  the  transmission  of  gamna 
rays  through  rectangular  ducts  from  basic  scattering  principles,  using 
largely  Che  "albedo"  concept.  In  so  doing,  there  are  two  ends  In  view. 
The  first  Is  to  lay  down  a  set  of  tentative  analytical  techniques  to  be 
used  by  nuclear  engineering  designers  concerned  with  this  problem, 
pending  further  refinement  of  the  art;  the  second  Is  to  provide  some 
theory  which  experimental  programs  carried  out  by  NCEL  (and  other 
organizations)  can  undertake  to  check. 
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It  may  also  be  possible  to  obtain  some  sound  knowledge  of  how  the 
radiation  attenuation  factors  vary  with  the  geometric  scale  of  the 
physical  layout.  Such  a  knowledge  may  be  useful  in  permitting  analysis 
of  large-scale  entranceways  through  small-scale  experiments.  Therefore! 
conmencs  on  the  scaling  relationships  will  be  provided  where  something 
meaningful  can  be  said. 


II.  AKALYSIS  OF  THE  BASIC  PROBLEM 
A.  The  Albedo  Concept 

The  basic  application  of  the  albedo  concept  to  this  sort  of  problem 
is  not  new  (see*,  for  example,  Ref.  1).  For  the  sake  of  the  reader  not 
too  familiar  with  Its  use  In  such  a  case,  the  following  review  Is  pre¬ 
sented. 

Consider  an  Ideallzec  basic  situation  , In  which  an  Isotropic  point 
source  emits  radiation,  and  a  detector  Is  situated  such  that  It  re¬ 
sponds  only  to  radiation  scattered  from  a  small  plane  area,  A.  By 
"small"  area,  we  mean  that  Its  linear  dimensions  are  much  less  than  r^, 
the  distance  from  source  to  the  area,  and  r^,  the  distance  from  the 
area  to  the  detector.  (See  Figure  1.) 


A  set  of  cartesian  coordinate  axes  Is  established  with  the.  origin 
within  A,  and  In  such  a  way  that  the  Z-axls  la  normal  to  the  plane  of 
the  area  A,  and  the  point  source  is  In  the  X-Z  plane.  In  addition  to 
the  radial  distance  r.  ,  the  position  of  the  source  Is  Indicated  by  the 
angle  from  the  normal''^..  Obviously,  there  Is  no  reason  to  expect  the 
detector  to  be  In  the  X-Z  plane  necessarily,  and  thus,  In  addition  Co 
the  radial  distance  angles  Sj  and  e  are  required  to  estab¬ 

lish  Its  position. 

If  appropriate  means  are  provided  to  eliminate  any  passage  of  ra¬ 
diation  from  source  to  detector  except  by  scattering  off  A,  and  If  there 
Is  no  attenuation  by  media  In  the  path  from-source-to-area- to-decector, 
the  detector  reading  Is  readily  seen  to  be  (Ref.  1,  page  335): 


D  AS  cos?, 
o 


(1) 


where  Do  Is  a  detector  response  in  the  direct  beam,  at  one  unit  of  dis¬ 
tance  from  the  source;  and  1 B a(?j is  the  "differential  direction¬ 
al  albedo"  (hereinafter  referred  to  as  the  "differential  albedo"), 
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or  reflection  factor,  for  the  radiation.  In  addition  to  the  geometric 
angular  factors,  the  differential  albedo  depends  also  on  the  energy  of 
the  incident  radiations  and  the  nature  of  the  scattering  material  undir 
the  surface  of  A.  It  gives  Che  proportion  of  radiation  incident  on  A 
which  is  reflected  off  in  the  given  direction,  per  ateradian. 

If  D  and  Do  are  numbers  of  particles  (or  photone),  a  must  be  a 
"number  albedo";  if,  on  the  other  hand,  the  detector  measures  radiation 
energy  flux,  D  and  Do  represents  energy  and  3  is  an  "energy  albedo"; 
finally,  if  D  and  Do  indicate  "dose,"  3  is  a  "dose  albedo." 

There  is  available  at  present  considerably  more  information  on 
"total  albedo,"  the  proportion  of  radiation  scattered  in  totof  than  on 
the  "differential  albedo,"  which  is  that  scattered  into  specific  direc¬ 
tions.  If  information  is  only  available  on  total  albedo,  one  must  make 
some  assumptions  as  to  the  actual  distribution  with  emergent  angles  92 
and  a.  For  example,  if  an  isotropic  assumption  la  made.  Equation  1 
becomes: 
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°o^*T 
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where  a«aaj(d^)  is  the  total  albedo  factor  for  the  radiation  energy 
and  reflecting  material  considered,  and 


■ J  3(e^,e^,e)  sinSj  de^  da  (2b) 

upper  hemisphere 


If  a  number  of  scattering  areas  of  different  orientation  are 
involved,  the  detector  dose  becomes: 
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B.  Approach  and  Assiuptlons 

The  basic  problem  to  be  solved  is  depicted  in  Figure  2.  This 
shows  a  two-legged  duct,  with  legs  intersecting  perpendicularly,  having 
a  rectangular  cross  section.  The  height  of  the  rectangle,  H,  is  the 
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dimension  perpendlculsr  to  the  plsne  of  the  paper  in  the  figure;  It  Is 
taken  to  be  the  same  for  both  legs  —  a  usual  situation.  The  trans¬ 
mission  of  radiation  through  the  duct  will  be  calculated  f*lrst  with 'the 
assumption  that  only  those  areas  which  can  be  seen  by  both  the  source 
and  the  detector  are  important.  The  ease  is  considered  to  be  the 
"basic"  ease;  its  solution  is  the  first  approximation  to  the  total 
solution;  the  scattering  areas  involved  are  called  "prime"  scattering 
areas  and  are  shown  as  areas  through  A^,  inclusive,  in  Figure  2, 
Areas  1  and  2  are  wall  areas;  areas  3  and  4  each  include  both  roof  and 
floor. 

Although  the  jlscusaion  is  made  with  gamna  rays  primarily  in  mind, 
the  formulas  derived  are  applicable  to  any  type  of  radiation  which 
reasonably  confonss  to  the  assumptions. 

At  this  basic  stage  in  the  development  of  the  theory,  scattering 
off  one  surface  only  will  be  considered.  Multiple  scattering  effects, 
to  the  degree  they  are  significant,  are  discussed  In  Section  V-C. 

A  certain  degree  of  approximation  will  be  necessary  in  order  to 
provide  some  mathematical  simplicity.  Correction  terms  on  the  order 
of  10  percent  or  greater  will  be  retained,  those  of  the  order  of  a  few 
percent  or  less  will  be  Ignored. 

Some  simplifications  are  possible  If  certain  dimensional  symme¬ 
tries  or  equalities  are  asausied  or  If  simple  expressions  for  the  albedo 
factors  are  assumed.  For  the  present  analytical  development  these  are 
not  assumed,  so  that  the  expressions  derived  will  be  generally  useful. 

It  is  assumed  that  the  leg  lengths  are  appreciably'  greater  than 
the  widths  and  the  height.  More  precisely,  the  length  of  the  shorter 
leg  should  be  st  lesst  three  or  four  times  the  width  of  the  wider  leg 
(or  height)  for  assured  sccurecy.  Also,  to  minimize  direct  trensmie- 
slon  of  radiation  from  source  to  detector,  both  legs  should  not  be  less 
than  s  few  feet  in  length. 

It  Is  assumed  that  the  leg  dimensions  are  much  less  than  a  mean 
free  path  of  the  radiation  In  air,  so  that  interaction  with  the  air 
may  be  considered  negligible. 

All  the  linear  dimensions  of  the  duct,  such  as  width  and  height, 
are  on  the  order  of  or  greater  than  a  mean  free  path  of  the  radiation 
In  the  material  of  which  the  walls  are  composed.  This  requirement,  for 
example,  would  I'ndlcate  that  the  smallest  duct  through  concrete  to 
which  the  analysis  Is  applicable  would  have  s  width  and  height  of  sev¬ 
eral  inches. 
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The  wall  naeerlala  are  unlfoni  In  eompoaltlon  and  dpnsltyt  eapa- 
dally  ar.  the  corner.  If  the  wall  Is  made  of  a  liner  of  some  material 
Imbedded  In  another  material,  the  liner  has  a  thickness  on  the  order  of 
or  greater  than  a  mean  free  path  of  the  radiation  in  the  liner  material. 

The  radiation  source  is  located  at  the  center  of  the  duct  entrance; 
the  detector  Is  located  at  the  center  of  the  duct  exit.  Under  the 
assumptions  given,  a  slight  variation  away  from  the  center  line  of  the 
duct  for  either  source  or  detector  should  not  cause  too  great  a  varia¬ 
tion  in  results.  Furthermore,  a  source  distributed  over  the  awuth  of  the 
duct  entrance  may  be  eonaidered  as  concentrated  in  the  assumed  spot 
without  severe  Inadcuracy,  provided  D  is  appropriately  calculated  for 
such  a  case,  ^ 

Most  of  the  above  restrictive  asauBq>tions  may  be  relaxed  to  some 
degree  without  introducing  serious  inaccuracies.  However,  the  exact 
degree  of  approximation  in  so  doiftg  is  not  thoroughly  studied  and  is 
largely  left  for  the  future. 

C.  Basic  Calculation 

The  basic  calculation  will  consider  scattering  from  areas 
through  A^,  inclusive.  (Areas  in  the  plane  of  the  paper  are  doubled  so 
as  to  account  for  both  top  and  bottom  surfaces.)  These  areas  are 
readily  obtained  from  geoewtrlcal  considerations  and  are  listed  in 
Table  I. 


The  values  of  eoss^,  as  requited  in  Equation  1,  are  also  readily 
obtained  and  are  given  in  Table  I.  In  oiaking  this  computation  for  each 
area,  the  radiation  is  eonaidered  to  strike,  on  the  average,  a  point 
whose  location  la  conveniently  expressed  smthematleally  and  which 
roughly  approximates  the  centroid  of  the  area. 

The  values  of  r,  and  T2  in  Equation  1,  as  well  as  the  symbol  de¬ 
noting  Che  proper  value  of  the  albedo  factor,  for  the  areas  in  question 
are  given  In  Table  II. 

From  these  quantities  and  the  use  of  Equation  1,  Che  doses  result¬ 
ing  from  the  various  areas  are  readily  arrived  at,  and  are  listed  in 
Table  III. 
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Table  II.  Distance  and  Albedo  parameters  In  Analysis  of 
Basic  Problem  and  Lip  Transmlaslon  Effect 


■Mjll 

Mean  Distance 

Albedo  Factor 

From  Source 

From  Detector 

H 

U  (1  +  Pj) 

‘1 

2 

Lj  (1  +  Pj) 

^2  ’ 

‘2 

3 

h 

■3 

4 

4 
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■4  (-  *3)  . 

5 

Li  (1  -  Pj) 

Lj  (1  +  Pj) 

•5 

6 

h  ®2> 

1.2  Cl  ■  P;) 

*6 

7 

4  (1  -  fJ2> 

8 

h 

Lj  (1  -  Pi) 

*8 

Note:  =  2  L2  ’  ^2  “  2  ’  ^3  "  2  L2 


ij 
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^bealc 


or 

^bea  Ic 


where 


a,, 
areas  1, 


total  doaa  expected  from  the  prime  scattering  areas  Is  summed 
contributions  of  the  four  areas  and  Is  given  as  follows; 


°o  Wj  H  jr  a^  Sj, 

2  ^  3  ^  2  I  1  +  +  ^2^ 


(3«) 


°i  h  h  h'> 


^1  -71 

^1 


h  -  rq 
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n  H 

2  L- 


<3b) 


y  (3c) 


_ t2 _ +  + 

^2  ^  ”  ^1  ^  ”  ^2 
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a2i  a^,  and  a^  are  the  differential  dose  albedos  from  prime 
2,  3,  and  4,  respectively. 
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It  li  to  b«  noted  ehet  0^  li  the  theoretical  doae  expected  (with¬ 
out  wall  scattering  effects)  at  the  comer,  that  Is,  where  the  mld-llnes 
of  the  two  legs  intersect.  It  Is  useful  to  separate  out  this  particular 
factor,  since  certain  effects  can  be  approxlSMted  by  appropriate  modifi¬ 
cations  to  this  term.  This  will  be  seen  later. 

It  la  obvious  from  a  glance  at  Equation  3b  that  the  last  two  factors 
are  nondlsMnalonal ,  and  that  for  a  given  geometrical  configuration,  their 
product  la  Independent  of  any  scale  selection.  Under  circumstances  which 
allow  this  basic  analysis  to  be  an  adequate  approximation,  the  sqallng- 
Independence  of  this  part  of  the  expression  permits  ready  adaptation  of 
experlisental  solutions  on  a  small  sMdel  scale  for  solution  of  large- 
scale  problems  of  this  nature.  (Possibly  Important  correction  factors 
may  make  this  more  difficult,  however,  as  we  shall  see  presently.) 


III.  CALCULATIONS  OF  COKBECTIOMS  FOR  VARIOUS  CORNER  LIP  EFFECTS 

A.  General 

Certain  important  other  effects  are  desirable  to  consider,  for  In 
some  cases  they  are  quite  lsq>ortant  compared  to  results  of  the  basic 
computation,  even  under  circumstances  In  which  the  assumptions  listed 
In  Section  II. B  are  smt.  Theae  effects  are  related  to  the  existence  of 
the  so-called  comer  "lip,”  and  chair  Inclusion  In  the  analysis  consti¬ 
tutes  a  sort  of  "first-order”  correction  to  the  basic  solution. 

B.  Corner  Lip  Transmission  Effects 

In  the  basic  calculation  it  Is  assumed  that  the  corner  lip  (the 
inner  edge  of  the  Intersection  of  the  two  legs)  Is  completely  opaque  to 
the  radiation.  This  can  never  be  precisely  true,  and  In  some  cases 
radiation  penetration  of  the  lip  can* be  quite  significant. 

If  it  is  assumed  chat  radiation  absorption  of  a  ray  passing  through 
the  lip  Is  exponential  in  character,  one  can  show  on  an  elementary  basis 
Chat  Che  amount  of  radiation  passing  through  the  lip  Is  the  same  as  If 
all  the  radiation  within  a  certain  cut-off  point  were  transmitted  and 
all  beyond  it  were  absorbed.  Furthermore,  the  distance  of  travel  within 
the  lip  material  of  that  particular  ray  passing  exactly  through  the  cut¬ 
off  point  can  be  shown  to  be  the  reciprocal  of  the  effective  attenuation 
coefficient  for  the  radiation  passing  through  it  (assuming  the  use  of  an 
"effective  attenuation  coefficient"  for  dose  absorption  Is  valid).* 
Appendix  I  shows  the  proof  of  these  statements. 


Ae  tlmplftst  assumption  for  an  affective  absorption  coefficient  Is 
that  which  accounts  for  all  energy-absorblng  processes  aiid  assumes  that 
all  energy  scattered  but  not  absorbed  Is  scattered  through  a  negligibly 
small  angle.  The  use  of  the  "energy  absorption  coefficient"  gives  thus 
a  reasonable  value  to  use  In  this  regard. 

On  this  basis,  then,  the  effect  of  comer-lip  penetration  can  be 
readily  approximated  by  an  Increase  In  the  scattering  area  beyond  the 
primary  scattering  areas.  The  areas  on  Figure  2  which  are  designated 
A5  through  Ag,  Inclusive,  show  the  new  scattering  surfaces  which  con¬ 
tribute  to  detector  response,  through  the  albedo  process.  The  calcula¬ 
tion  of  the  values  of  these  areas  and  other  variables  required  to  deter¬ 
mine  their  effect,  according  to  Equation  1,  is  straightforward  and  is 
susnarlsed  In  Tables  I,  II,  and  III,  along  with  previously  mentioned 
data  for  the  basic  area  contributions. 

It  Is  convenient  to  break  up  this  lip  transmission  effect  contri¬ 
bution  Into  two  parts!  one  in  which  wall  scattering  occura  before  lip 
transmission,  and  the  other  In  which  lip  transmission  occurs  before  vail 
scattering.  We  find  then  that 
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^  O  ..  ^  .  n  A  J 


2  2  Lj  (I  -  3^)^  (1  <•  $2)’ 


(4a) 


°1  h  ®tl 


(4b) 


'tl 


(1  -  Pj^)  Sj  +  2  a^ 

2  L2  (1  -  (1  -  02>^ 


(4c) 
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and 


®o  *2  "  f  "  ^2^  *6  ^  ^2  *8 


'tr2 


*1  "2  "  f 
,  3  ,  2  1 


2  L^-  Ij*  "  2  02  <1  “  Pi)^  <1  -  ^2^^ 


<5n) 


®1  ^1  ^2  ^3^  °t2 


(5b) 


't2 


<1  -  0^)  «6  ^  ^  *8 
2  02  (1  -  0^)^  (1.-  02>^ 


<5c) 


Dj^  has  baan  pravioualy  daflnad  by  tha  sat  of  Equations  3e.  Is  tha 
enargy  absorption  eoaffieiant  of  tha  priaary  radiation,  is  tha 
coaffieiant  of  tha  radiation  raflsctad  from  ths  surfaeas.  Sae  Table  IV 
for  nuaarieal  data  on  thasa  eoafficiants.  * 

It  it  to  ba  noted  that  in  Equations  4b  and  5b  tha  two  latter  terms 
of  each  aquation  ara  nond loans lonal;  however,  the  product  of  these  two 
terms  is  not  Independent  of  tha  gaomatrlc  scale  selection  but  Is  in¬ 
versely  proportional  to  tha  scale  siae.  We  dote  that  the  attenuation 
coefficients  are  fixed  and  do  not  change  value  with  a  change  in  geometric 
scale.  Since  tha  scaling  relationships  for  similar  factors  in  the  basic 
problem  and  this  particular  correction  are  not  the  same,  the  problem  of 
scaling  up  from  model  experiments  to  larger-scale  prototypes  is  seen  to 
be  impossible  unless  some  way  is  found  to  separate  experimental  data 
into  the  basic  part  and  the  lip-effect  terms,  of  unless  the  lip-effect 
terms  can  be  shown  to  be  negligible  at  all  scales  considered. 

Djyj  and  added  to  give 


D 


tr 


(4  0j  02  03) 


(6a) 


where  °t  ■  ®tl  ®t2 


<6b) 
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Table  IV.  Nuclear  Quantities  Used  In  Dflct  Attenuation  Formulas 


S  (Mev) 

0 

Ua 

• 

H,  (cm~S 

/  Z  N  N  2 

{1:7 

Material 

0.50 

.0473 

.0473 

17.00 

1.25 

.0473 

.0434 

20.10 

Earth 
(100  pcf) 

6.00 

.0473 

.0300 

45.70 

0.50 

.Q695 

.0695 

11.55 

1.25 

.0695 

.0630 

13.70 

Concrete 
(145  pcf) 

6.00 

.0695 

.0442 

31.30 

0.50 

.3140 

.2240 

3.38 

1.25 

.2450 

.1900 

4.68 

Iron 

(475  pcf) 

6.00 

.2260 

.1742 

5.47 

0.50 

6.260 

1.120 

1.68 

1.25 

2 . 545 

0.375 

1.51 

Lead 

(705  pcf) 

6.00 

1.318 

0.422 

1.19 

Note  #1.  The  absorption  coefficients  ()i)  for  earth  and  concrete  are 
relatively  constant  in  the  ranges  0.2  -  0.6  Mev  and  4-10  Kev,  which 
cover  the  significant  ranges  of  energy^  for  single  90-degree  scattered 
and  initial  bomb  radiation  spectra,  respectively. 

Note  #2.  The  absorption  coefficients  are  Inversely  proportional  to  the 
density  of  toaterlal  with  roughly  the  same  average  atomic  number.  For 
densities  different  than  those  listed  use  the  ratio;  up/pg,  where  p  Is 
the  actual  density  and  pg  the  density  listed  In  this  table.  The  same 

thing  Is  true  of  the  term  ( i  r 

V  ® 
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C.  Corner  Lip  In-Scatterlng  Effect 


Not  only  may  the  corner  lip  transmit  some  of  the  radiation  photons 
(or  particles)  eosipletsly,  hut  It  may  also  serve  to  scatter  sosM  of  them 
one  or  more  times  In  the  passage.  Such  scattering  will  In  part  redirect 
radiation  toward  the  detector.  Thus  the  detector  not  only  "sees'*  the 
prime  and  additional  scattering  areas  on  the  outer  walls  by  means  of 
their  radiation  scatter,  but  also  "sees"  the  comer  lip  as  a  "bright" 
source  —  almost  a  line  source. 

In  analyzing  the  contribution  of  this  effect  we  make  certain  sits* 
pllfleatlons.  In  similar  spirit  to  that  used  in  the  preceding  section 
we  utilise  an  "effective  attenuation  coefficient"  for  radiation  passing 
to  and  from  each  scattering  center,  thus  again  using  what  might  be 
considered  a  "straight-ahead  approximation."  The  previous  approach  la 
modified  to  the  extent  that  we  recognise  a  small  but  definite  probabil¬ 
ity  of  scattering  Into  the  direction  which  will  cause  the  radiation  to 
hit  the  detector.  Thus  we  have  a  single-scatter  approximation  for  the 
most  part,  but  the  approach  docs  not  eliminate  multiple  scattering  pro¬ 
vided  all  the  scattering  processes  but  one  are  considered  of  a  small 
angle  niture.  Scattering  through  two  or  more  appreciable  angles  of 
large  amount  is  considered  of  negligible  proportion  (an  approximation 
which  is  admittedly  not  entirely  valid  for  photons  or  particles  of  quite 
low  energy). 

The  computation  of  the  effect  for  gamma  ray  photons  Is  based  on  the 
Kleln-Nlahlna  scattering  fonsula  (Ref.  2).  Its  use  In  this  particular 
case  Is  explained  in  detail  in  Appendix  II.  The  results  of  the  detailed 
analysis  are  as  follows: 


W,  W,  H  Z  N  K 

I  .  - : — ; — ::: -  (7a) 

3  3  2  3  3 

4  l/  (1  -  (1  - 

_  Dj  (4  $2  ^3^  °s 

,  Z  N  K 

8  *  z  - - r 

2  L2  (1  -  Pi)-*  (1  -  P2) 


where 


(7c) 


And  whAtA  Z  m  tiumbar  of  eloetrons  per  atom  of  Che  seaccerlng  material 

M  a  number  of  aeons  per  unit  volume  of  the  acaCterlng 
material 

lit  ■  K(d  •£  )  Is  the  Klein- Nlshlna  coefficient  for  scattering 
probabfilty  per  electron  (see  Apprendlx  II  for  detailed 
definition  of  this  term  and  Its  arguments  -  also  note 
Figure  3  and  Table  IV) 

fl,  ■  90®  -  Oj^  -  0^ 


2 

From  Figure  3  one  can  obtain  the  value  of  K/ro  for  various  values 
of  dg  and  So*  From  Table  IV  one  can  obtain  appropriate  values  of 
Z  N  The  multiplication  of  corresponding  values  of  these  two 

factors  Will  give  values  for  Z  N  (The  fact;or  drops  out  of 

the  final  expreaelon  as  superfluous;  It  Is  Included  herein  for  reasons 
related  to  the  basic  derivation.  The  meaning  of  rp  la  "the  classical 
radius  of  the  electron";  its  value  Is  about  2,8  x  ItT^^  cm.)  The 
geometrical  significance  of  and  0^  Is  Indicated  In  Figure  IS  of 
Appendix  II. 

We  see  from  Equations  7  that  dimensional  scaling  has  the  same  effect 
on  this  cpntrlbutlon  as  It  has  on  the  contribution  of  the  corner  lip 
transmission  effect,  so  all  the  corner  lip  effects  scale  In  similar 
fashion,  which,  as  we  have  pointed  out,  Is  not  the  same  as  for  the  basic 
effect  of  the  prime  scattering  surfaces. 

It  also  la  to  be  noted  that  the  contribution  of  this  effect  Is 
roughly  Inversely  proportional  to  the  value  of  and  therefore 
decreases  with  Increasing  atomic  numbers  In  much  the  same  way  as  In  the 
comer  lip  transmission  effect.  (Since  Is  roughly  proportional  to 
NZ,  one  of  the  two  factors  in  the  term  In  the  denominator  Is  can¬ 
celled  in  the  proportionality  calculation.) 
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Figure  3.  Klein^Nlshlna  coefficient  for  scattering  probability  per  electron. 


IV.  SUMHARY  OF  METHOD  OF  ANALYTICAL  CALCULATIONS 


A.  Attenuation  Factors 

'  He  are  now  at  the  point  of  being  able  to  sunmarlse  the  analytical 
formulas  derived  up  to  this  point,  with  the  hope  that  answers  obtained 
for  common  situations  meeting  the  qualifications  given  in  Section  II3 
will  be  reasonably  correct,  at  least  for  ganaa  rays. 

Adding  Equations  3,  6a,  and  7b,  we  get  the  total  dose  at  the  end  of 
the  second  leg: 


°TOT  *^1  ^1  ^2  ^3^  ‘^TOT 


(Sa) 


where 


Sot  -  °b  +  °t  +  S 


(8b) 


We  'are  now  at  the  point  where  we  can  Introduce  the  attenuation 
factor  and  solve  for  it.  For  the  pturpose  of  this  paper,  we  will  use  the 
following  definitions; 


F-  se  D_/D  ,  being  Che  ratio  of  dose  at  the  end  of  the  second  leg  to 
°  chat  at  the  reference  point,  one  foot  from  the  source 
in  a  nonabsorptlve,  nonscattering  medium 

F.  as  D./D^,  being  Che  ratio  of  the  dose  at  the  intersection  of  the 
^  leg  center  lines  to*thaC  at  the  reference  point,  one 
foot  from  Che  source  in  a  nonabaoFptive ,  nonscattering 
medium 

F,  ■■  P../O1  >  being  the  rstlo  of  the  dose  at  the  end  of  the  second 

leg  to  that  at  the  intersection  of  the  leg  center  lines 


He  will  commonly  call  F]^  the  attenuation  factor  for  the  first  leg,  and 
Fj  the  attenuation  factor  for  the  second  leg.  Obviously,  F^  •• 


He  now  find  it  possible  to  consider  the  attenuation  factors  for  the 
two  legs  separately.  We  readily  see  that  ; 


Fi  . 


(9) 


and  that 


Fj  «  4  0^  02  ^3  Sot 


(10) 
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It  la  polntad  out  that  in  our  dlacuaalon  of  modal  icallng  ralatlon- 
ships  in  Saetlons  II. C,  III.B,  and  III.C,  the  discussion  was  relative  to 
the  second-leg  attenuation  only,  since  the  factors  involving  first-leg 
attenuation  were  eliminated  from  consideration. 

By  separating  our  factors  in  this  way,  we  also  allow  ourselves,  if 
we  wish,  to  use  a  different  reference  point  for  overall  attenuation 
(rather  than  the  point  at  which  Dq  Is  measured)  without  affecting  the 
analysis  for  which  is  the  store  difficult  part. 

B.  Format  and  Bata  for  Calculations 

I 

Tables  V(a)  and  V(b}  are  a  fonut  which  will  permit  ready  calcula¬ 
tion  for  this  type  of  problem,  using  the  information  and  formulas  pro¬ 
vided.  Since  F^  is  readily  calculable,  the  forms  are  designed  for  the 
purpose  of  calculating  F^. 

Since,  at  the  present  state  of  the  art,  precise  knowledge  of  direc¬ 
tional  albedo  la  usually  lacking,  whereas,  on  the  other  hand,  total 
albedo  Information  is  more  generally  available  (Refs.  3,  4,  and  5), 
some  sort  of  assumption  about  directional  distribution  must  be  made. 

The  fonss  arc  drawn  up  to  permit  use  of  the  Isotropic  assumption.  It  Is 
readily  determined  (see  Equation  2b)  that  under  Che  isotropic  asstmiptlen 
the  differential  albedo  equals  the  total  albedo  divided  by  2n.  Values 
for  concrete  and  low-energy  photons  are  given  In  Figure  4. 

Table  IV  gives  other  various  tabulated  information  needed  in  the 
solution  of  this  problem,  for  various  given  energies.  The  energies 
used  as  arguments  are  selected  for  the  following  reasons!  6.0  Mev 
represents  the  energy  of  the  most  penetrating  (an4  frequently  the  most 
abundant)  photons  In  the  Initial  radiation  from  nuclear  weapons  • 
explosions  in  air.  1.2S  Mev  Is  the  average  energy  for  eobalt-60  gamma 
radiation,  often  used  for  experimental  work;  it  la  also  a  convenient 
energy  to  use  Co  analyze  fallout  radiation  problems.  0.3  Mev  is  close 
to  the  value  of  photon  energies  after  a  90-degree  Compton  scattering, 
for  Incident  ganma  energies  of  general  interest,’  alSo  it  Is  the 
approximate  energy  of  poaltron  annihilation  radiation,  which  contrlbuCea 
heavily  Co  the  albedo  from  6-Mev  Incident  rays. 

It  might  be  pointed  out  that  the  data  are  not  highly  sensitive 
over  wide  regions  to  precise  photon  energies.  For  example,  the  energy 
absorption  coefficient  for  most  materials  la  fairly  constant  from  ^  to 
10  Mev,  and  thus  the  data  for  6  Mev  can  be  used  for  such  cases.  This 
coefficient  Is  also  quits  constant  for  concrete  and  earth  from  0.2  to 
0.6  Mev,  hence  the  use  of  0,5  Mev  as  representative  of  scattered 
radiation  energy  Is  sufficiently  accurate. 
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Table  V(a).  form  for  Dose  Computation 


“  1  +  Pj  (1  +  Pp  1 .  p^  I  -  Pj 

-(  )+(  )+(  )+(  )  G^» 

'  0 


*  Note;  Use  of  Table  IV  requires  all  lengths  to  be  measured  In  centi¬ 
meters,  all  values  of  linear  absorption  coefficient  to  be  In  reciprocal 
centimeters . 


2 


Albedo 


Figure  4;  Differential, directional  doae  bedo  for  1.25  and  0.5  Hev  gainna 
photons  froa  concrete  (Isotropic  asaunptlon) . 
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Tabls  VII.  Dl£f«rentl«l  DlrectlOMl  Doae  Albado  for  Gaana  Kays 
Scattcrsd  Fron  Gonerata  (Xaotropie  Aaaumpelon) 


Coa6^ 

E  «  0.5  (Mcv) 

.  0 

E  -  1.25  (Mav) 

0 

0.01 

.0628 

.0522 

0.02 

.0588 

.  .0509 

0.03 

.0572 

.0493 

.  0.04 

.0552’ 

.0469 

O.OS 

.0547 

.0453 

0.06 

.0530 

.0437 

4 

0.07 

.0518 

.0424 

0.08 

.0517 

..0413 

0.09 

.0495 

.o4o2 

0.10 

‘  .0491 

.0390 

0.20 

.0394 

.0287 

0.30 

■  .0326 

.0218 

0.40 

.0275 

.0186 

0.50 

.0232 

.0130 

0.60 

.0189 

.0105 

0.70  . 

.0164 

a 

.00843 

0.80 

.0146 

a 

.00684 

0.90 

.0132 

.00588 

1.00 

.0118 

.00508 

V.  TECHNiqUES  FOR  HANDLING  CERTAIN  COMPLICATIONS  AND  VARIATIONS 

A.  General  Comment 

Tl^ere  are  certain  conslderatlona  In  the  above  formulation  which 
have  been  alighted,  for  the  sake  of  avoiding  too  great  a  degree  of 
elaboration  and  because  adequate  techniques  for  handling  them  to  the 
desired  degree  of  accuracy  are  not  readily  available  at  present.  Under 
certain  clreumstancea ,  however,  they  may  become  somewhat  significant, 
providing  correction  factors  In  excess  of  10  percent.  Likewise,  If  we 
encounter  problems  In  which  the  given  conditions  and  physical  parameters 
do  not  meet. those  specified  In  Section  II. B,  certain  changes  In  the 
approaches  niay  be  desirable. 

Even  though  some  of  these  modifications  may  be  difficult  to  handle 
analytically  to  the  desired  precision,  we  can  give  some  recommendations 
as  to  how  they  can  be  appropriately  accomplished,  at  least  approximately. 

B.  Variation  In  Source  Configuration 

In  practical  cases  it  Is  unlikely  that  one  will  be  dealing  with  a 
source  concentrated  at  the  center  of  the  entrance  to  the  duct  system. 

One  la  much  more  likely  to  find  the  source  distributed  over  the  entrance 
or  to  find  the  radiation  streaming  In  from  outside.  If  such  radiation 
passes  down  the  first  leg  In  an  approximately  symmetrical  way,  one  might 
expect  the  foregoing  analysis  to  be  generally  valid  as  regards  the 
scattering  from  the  comer  walls  and  the  comer  lip  effects.  We  would 
expect  the  previously  derived  formulas  for  F2  reasonably  valid 

in  such  ease. 

On  the  other  hand,  the  value  for  F^  may  be  rather  strongly  depend¬ 
ent  upon  the  character  of  the  source  at  the  entrance,  for  two  reasons: 

The  angular  distribution  of  the  Incoming  radiation  may  not  be  such  as 
to  lead  to  a  simple  Inverse-square  relationship  of  dose  rate  with 
distance;  and  a  more  appropriate  reference  point  may  be  necessary  or 
desirable  as  the  basis  for  the  attenuation  factor  down  the  first  leg. 

The  question  arises,  then,  for  each  source. situation;  What  Is  the 
value  of  the  dose  rate  at  a  dlstrance  along  a  rectangular  tunnel,  or 
hole,  as  compared  with  the  dose  rate  at  the  entrance  or  at  some  other 
appropriate  reference  point?  The  following  subsections  discuss  some 
of  the  more  Important  of  these  situations. 
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1,  Ptstrlbutad  aourca  at  duct  entrance 

In  this  type  of  source  we  assume  that  the  radiation  source  Is  at 
the  entrance  of  the  duet,  hut  Is  equally  distributed  across  It.  This  Is 
what  is  normally  called  the  "Isotropic"  case  (Ref.  6,  Sect.  4.12).  In 
this  case  the  readings  down  the  first  leg  are  proportional  to 
ln(l  +  W^H/k  x^),  provided  the  width  and  height  of  the  first  leg  do 
not  differ  too  radically  in  amount.  Then  the  attenuation  factor 
depends  on  what  Is  used  aa  the  reference  point.  If  the  reference  point 
Is  taken  aa  one  unit  of  distance  away  from  the  center  of  the  entrance 
face' of  the  duct,  the  first  leg  attenuation  factor  becomes,  approxi¬ 
mately) 


2  2 

since  ln(l  +  W.B/xL.  )  approaches  W. R/xl.,  for  an  argument  much  greater 
than  unity. 

If  the  reference  point  la  taken  as  being  three  feet  above  a  smooth, 
infinite  plane,  covered  to  the  same  amount  per  unit  area  with  distrib¬ 
uted  radiation  sources,  emitting  isotropically,  one  can  show  that  (see 
Appendix  III): 


..  W,  H 

V  at _ i _ 

*^1  2 
'•  34  1.^ 


2.  Cosine  emission  source 

If  the  radiation  crossing  the  plane  of  the  duet  entrance  has  an 
angular  distribution  such  that  the  amount  passing  through  the  entrance 
from  directions  at  various  angles  to  the  leg  axis  Is  proportional  to  the 
cosine  of  the  direction  line  with  the  leg'axls,  the  distribution  Is 
called  "cosine."  It  Is  Interesting  to  note  that  this  situation  can 
actually  be  created  by  a  superposition  of  parallel,  broad- beam  radiation 
beams  coming  with  equal  Intensity  from  each  direction  In  toward  the  duct 
opening.  It  Is  thus,  with  respect  to  an  Isotropic  Instrument  outside 
the  duct  entrance.  Isotropic;  but  as  far  as  the  amount  penetrating  the 
duct  entrance.  It  Is  angle- dependent,  since  the  amount  penetrating 
depends  on  the  angle  of  "aspect"  of  the  entrance  area  to  the  radiation 
coming  from  a  certain  direction. 
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Under  these  elrcumdtences,  it  may  be  shown  thet  the  ettenuetion 
factor  In  the  first  leg,  based  upon  a  reference  value  at  or  Just  outside 
the  entrance  to  the  duct.  Is  <Re£.  6,  par.  4.12): 


i _  (13) 

2 

2it  L, 

,  *1 

provided  the  ratio  between  the  first  leg  width  and  height  does  not  vary 
too  laerkedly  from  unity. 

Under  auiny  general  conditions,  a  crude  approximation  to  the  problem 
of  Initial  nuclear  weapon  radiation  entering  a  shelter  entranceway  may 
make  use  of  this  approach. 

3.  Othar  typea  of  sourcea 

The  nature  of  the  incoming  source  radiation,  particularly  its 
directional  distribution,  can  make  a  great  difference  In  the  value  of 
7]^.  ?or  example,  monodirectlonal  radiation  coming  In  parallel  to  the 
axis  of  the  first  lag  will  reach  a  detector  at  the  end  of  the  first  leg 
unattenuated.  In  aueh  a  case  the  value  of  Pi  Is  obviously  unity.  On 
tha  othar  hand,  a  monodirectlonal  bean  parallel  to  the  entrance  of  the 
duct  system  will  not  enter  at  all,  and  the  value  of  for  this  case 
la  sero. 

If  the  directional  distribution  of  radiation  coming  from  without  Is 
Irregular,  one  may  consider  the  proportion  of  the  radiation  coming  In 
through  a  solid  angle  subtended  by  the  entrance  to  the  duct,  as  seen 
from  the  point  of  Intersection  of  the  two  leg  center  lines,  compared  to 
Che  total  coming  from  all  directions  to  a  detector  on  the  outside.  This 
proportion  will  give  the  proper  value  of  based  on  the  said  outside 
reference  point.  The  situation  covered  In  Subsection  2,  above.  Is  a 
simple  case  for  this  situation.  On  the  other  hand,  If  Che  radiation 
Is  coming  predominantly  from  one  side,  this  approach  cannot  be  used 
blindly.  The  major  source  of  radiation  down  the  first  leg  In  such  a 
case  might  well  be  that  scattered  off  a  wall  Just  Inside  the  entrance 
to  Che  first  leg. 

C.  Multiple  Scattering  from  Duct  Leg  Walls 

Hitherto,  we  have  assumed  that  the  radiation  detected  at  the  end 
of  the  first  leg  is  that  coming  directly  from  or  through  the  entrance 
to  the  duct;  likewise  we  have  assumed  that  the  radiation  reaching  the 
detector  at  the  end  of  the  second  leg  Is  that  with  has  been  contributed 
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only  by  reflection  from  aurfeeee  through  Ag,  Incluelye,  plus  that 
iiweeettered  by  the  corner  lip.  It  Is  obvious,  however,  that  such  an 
asstaptlon  Ignores  the  possibility  of  an  appreciable  contribution  from 
duct  leg  walls  outside  those  already  designated  as  prlauiry  or  llp- 
tranamlaslon  scattering  surfaces.  The  question  arises:  What  error  Is 
Involved  In  this  assumption? 

It  smst  be  confessed  that  under  many  circumstances  within  the 
asaueptlona  of  this  paper  there  are  probably  situations  In  which  the 
error  Involved  In  Ignoring  multiple  scatter  la  greater  than  10  percent. 

Unfortunately,  this  subject  has  not  been  too  well  analysed  In  the 
past,  emcept  for  conditions  In  which  the  albedo  scattering  Is  independ¬ 
ent  of  the  Incident  angle  and  depends  In  very  slaiple  (Isotropic  or 
cosine)  fashion  upon  the  angle  of  emergence.  Such  an  approach  ouiy  be 
better  for  neutrons  than  gasBc  fays  (see  Refs.  1,  3,  and  6),  and 
furthermore  It  Is  dependent  upon  having  good  albedo  Information.  Rough 
calculations  and  experimental  data  Indicate  that  this  affect  may  In¬ 
crease  leg  attenuation  factors  by  well  over  10  percent  for  gamma  rays 
passing  through  concrete  duets  (Ref.  3).  (Note  that  by  an  Increased 
attenuation  factor  we  are  Indicating  a  poorer  degree  of  attenuation.) 
Attenuation  In  the  second  leg  may  be  more  a^laeted  than  In  the  first 
leg,  since  the  albedo  for  light  materials  such  as  concrete  Increases  as 
the  energy  is  lowered  (Refs.  4  and  S). 

D.  Direct  Transmission  Effect 

For  ducts  of  smaller  slse.  It  Is  possible  that  appreciable  radia¬ 
tion  auy  be  transmitted  directly  from  source  to  detector.  The  compu¬ 
tation  of  this  effect  is  quite  straightforward.  We  cannot  apply  the 
results  directly  to  the  various  attenuation  factors,  .but  must  add  doses 
or  dose  rates  (absolute  or  relative). 

From  Figure  2  we  readily  see  that  the  distance  between  the  source 
and  detector  Is  -f  L7^.  However,  the  direct-line  thickness  of  wall 

material  is  (1  -  -  Pj)  CWe  assume  the  intervening 

material  to  be  homogeneous  —  if  not,  appropriate  modification  to  the 
analysis  can  be  made.) 


IC  appears  suitable  to  use  the  build-up  factor  as  determined  for  an 
ideal  case  in  which  source  and  detector  are  imbedded  in  Infinite,  homo¬ 
geneous  material  (see  Ref.  7).  We  thus  conclude  that; 
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where  ^  is  the  total  absorption  coefficient  for  gamaa  rays  of  the 
initial  energy  at  the  source  passing  through  the  given  duct  wall 
material.  The  build-up  factor,  B,  Is  a  function  of  the  number  of  mean 
free  paths  through  the  wall  material.  The  number  of  mean  free  paths  is 

equal  to  (1  -  -  ^2^V^1^  "**  ^2^  ’  obtained  from  Ref.  7, 

appropriate  extracts  from  which  are  provided  In  Table  VI. 


Table  71.  Dose  Build-up  Rectors  for  Aluminum  and 
Other  Light  Materials  (Rrom  Ref.  7) 


Photon 

No.  of 

Mean 

Free  Paths 

Energy 

(Mev) 

1 

2 

4 

• 

7 

10 

15 

20 

0.5 

2.37 

4.24 

9.47 

21.50 

38.90 

80.80 

141.00 

1.0 

2.02 

3.31 

6.57 

13.10 

21.20 

37.90 

58.50 

2.0 

1.75 

2.61 

4.62 

8.05 

11.90 

18.70 

26.30 

3.0 

1.64 

2.32 

3.78 

6.14 

8.65 

13.00 

17.70 

4,0 

1.53 

2.08 

3.22 

5.01 

6.88 

10.10 

13.40 

6.0 

1.42 

1.85 

2.70 

4.06 

5.49 

7.97 

10.40 

note:  The  source  Is  considered  point- Isotropic . 
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VI.  COHPARISON  WITH  EXFERIMEHT 


C<  Elsenhauer  has  perfonoed  some  experiments  with  small  concrete 
ducts  and  cobalt-60  (Ref.  8).  Figure  5  is  b  comparison  of  Elsenhauer 's 
experimental  results  and  the  LeDoux- Chilton  analysis  using  the  Isotropic 
differential  albedo  assumption.  Because  of  the  assumptions  made  In 
Section  XI. B,  correlation  with  experiments  of  small  L/W  ratios  was  not 
expected.  Overall  agreement  Is  good. 

In  order  to  measure  the  effect  of  the  corner  contribution, 
Elsenhauer  placed  a  lead  block  in  the  corner  of  the  duct  and  then 
meaaured  the  remaining  gamma  dose.  Figure  6  Is  a  camparlson  of  experi¬ 
mental  data  with  the  LeDoux-Chllton  analysis  for  this  case.  In  order 
to  make  the  calculations  for  this  case,  the  proper  values  for  lead 
from  Table  III  must  be  used  In  computing  Gt  and  Gs ,  which  are  the 
contributions  from  the  comer.  Elsenhauer  does  not  consider  his  data 
highly  accurate,  and  the  very  close  agreement  Is  considered  rather 
fortuitous  under  the  circumstances. 

Charles  Terrell  has  performed  experiments  on  large  (6- ft  by  6- ft) 
ducts  (Ref,  3).  Figure  7  Is  a  comparison  of  the  LeDoux-Chllton  analysis 
with  his  results.  Agreement  Is  again  good. 

For  these  analytical  comparisons,  none  of  the  complications  dis¬ 
cussed  In  Section  V  were  allowed  for. 

From  these  comparisons  It  thus  appears  that  a  basic  theoretical 
analysis  can  give  results  which  check  very  well  with  experimental 
results.  The  isotropic  distribution  assumption  for  differential 
directional  albedo  seems  to  be  an  excellent  choice.  The  effect  of 
using  more  accurate  albedo  data  awaits  the  publication  of  such  data. 
Meanwhile,  for  the  cases  considered,  the  LeDoux-Chllton  method  of 
analysis  appears  to  be  sufficiently  accurate  for  practical  shielding 
computations,  Data  from  Table  VXI  may  be  used  for  practical  shielding 
computations  based  on  the  isotropic  albedo. 


VII.  NEED  FOR  FUTURE  WORK 

The  Monte-Carlo  calculations  now  being  done  (such  as  the  work  In 
progress  by  D.  Raso  of  Technical  Operations,  Inc.)  should  yield 
sufficient  data  for  the  differential  dose  albedos  for  energies  from 
0.2  to  over  10  Mev,  for  the  materials  of  shielding  importance.  It  is 
hoped  that  this  work  will  permit  a  more  accurate  comparison  of  this 
analytical  approach  and  experimental  results. 


More  experimental  data  is  needed  for  distributed  sources,  for  other 
source  energies,  and  for  neutrons ,  Terrell,  of  Armour  Research  Founda- 
^  tlon,  Is  doing  work  In  this  direction.  Once  a  good  correlation  la 
obtained  It  will  be  poseihle  to  analyse  duct  and  entranceway  shelter 
radiation  streaming  problems  with  adequate  accuracy. 

There  Is  also  the  possibility  of  experimental  analyses  of  even  more 
complex  situations  through  work  with  ducts  of  small  scale. 
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Dose,  mr/hr 
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mx/hx 


Lj/W 


Figure  7«  Conpertson  of  the  LeOoux-Cbllton  enalysle  with  ASF 
experiment  using  cobalt- 60  and  a  concrete  duct, 

W  -  6  feet. 
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Appendix  I 

CORNER  LIP  TRANSMISSION  EFFECT 


The  case  to  be  considered  Is  a  rectangular  duet  as  shown  In 
Figure  8.  The  height  of  the  duct  Is  H,  the  dimension  perpendicular  to 
the  plane  of  the  paper.  It  Is  not  reasonable  to  assume  different  heights 
of  the  two  duct  legs,  since  the  simplest  way  of  Joining  the  duct  legs  Is 
to  require  a  comnon  height  H.  The  widths  of  the  ducts  are  and  W2. 

He  will  conalder  a  general  case  of  nonperpendlcular  duct  legs,  and  will 
particularise  to  the  perpendicular  case  later. 


Figure  9  llluatrates  Che  geometry  of  the  corner  lip  tranamlsslon 
effect  for  the  ray  which  Is  transmitted  through  the  lip  prior  to  being 
reflected  off  a  duct  wall.  For  convenience,  this  effect  wilt  hereinafter 
be  called  simply  the  "tranamiaslon  effect."  The  duct  has  legs  of  lengths 
and  L2«  Intersecting  at  an  angle  5.  We  consider  a  photon  from  the 
source  S  which  traverses  the  comer  through  a  thickness  before  being 
reflected  from  Leg  2  Into  the  detector. 


If  we  focus  our  attention  on  the  corner  lip  Itself,  Figure  10,  It 
can  be  seen  chat  there  Is  approximately  a  parallel  beam  of  photons  strik¬ 
ing  Che  comer  at  an  angle  Oh  (assuming  L]^»W]^).  If  Che  Intensity 
striking  Che  comer  Is  I(Bev/cm^) ,  and  the  thickness  of  the  beam  la  dx, 
then  Che  total  energy  per  unit  height  of  the  duct  would  be  1  dx.  Because 
of  absorption  within  the  material,  the  total  energy  per  unit  height  and 
beam  thicknes's  dx,  emerging  from  the  corner  would  be  I  e^MF  dx.  Since 
a  small  scattering  angle  and  multiple  scattering  may  In  part  still  pro¬ 
vide  photons  In  the  direction  of  Interest,  only  the  net. energy  absorp¬ 
tions  should  be  considered.  Hence,  the  ^  we  recotanend  to  be  used  here 
is  |i.,  the  linear  energy  absorption  coefficient,  closely  equivalent  to 
the  Effective"  absorption  coefficient. 


Since  only  the  gamma  rays  which  emerge  from  the  corner  lip  with 
no  (or  almost  no)  Interaction  are  of  Interest,  we  could  replace  the 
actual  situation  depicted  In  Figure  10  with  an  equivalent  transparent 
comer  as  shown  In  Figure  11.  This  "equivalent  lip"  Is  characterized 
by  a  dimension  T2  such  that  the  amount  of  photon  energy  which  Is  trans¬ 
mitted  by  this  truncated  lip  Is  the  same  as  that  transmitted  through  the 
actual  lip.  It  is  assumed.  In  the  equivalent  case,  that  any  photons 
which  strike  the  lip  material  so  as  to  to  pass  beyond  T2  are  totally 

absorbed. 
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gaoBBa  flux  I  through  corner 


Mathematically  this  can  be  expressed  by  the  following  equation; 


/ 


I  e 


-n,y 


dx 


(I-l) 


Trom  Figure  lOt 

y  ■  X  cot(B-aj)  +  X  cotOt^  (I"2) 

where  x  Is  the  distance  from  the  comer  to  the  path  of  the  photon. 
The  results  of  Integrating  Equation  I-l  are; 


X 


j^cotOj  +  cotCB-Oj^)^ 


(1-3) 


From  Figure  11; 


X  .  Tj  sln(S-aj) 


(1-4) 


Combining  Equations  1-3  and  1-4,  the  lip  dimension  Is; 

slnoL 

T, - —  (1-5) 

li  sinB 

fl 

In  a  slmllsr  fashion  the  dimension  of  the  lip  along  leg  1  for  the 
ray  which  strikes  the  wall  of  leg  1  first  and  then  penetrates  the  lip 
la: 

slna, 

T,  .  (1-6) 

U  8ln5 
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Fran  the  figures  it  otn  be  seen  thet  depends  on  Eq,  the  energy  of  the 
source,  end  depends  on  E^,  the  energy  sfter  reflection  from  the  wells 
of  the  first  leg  of  the  duet. 

Using  the  velues  of  end  the  various  velues  of  peremeters 
needed  for  computation  of  the  dose  from  these  transmission  areas  can  be 
coBiputed  and  are  listed  in  Tables  I  and  II  In  the  main  text.  It  should 
be  remembered  that  the  assumption  of  L»  W  still  holds.  This  Is  equi¬ 
valent  to  saying  that  sins  tanfit. 

If  6  Is  90  degrees,  then  Tj  and  T2  can  be  particularised  as  follows: 


(1-7) 


(1-8) 


slna, 

T.(90®)  . - 1 


TjOO**) 


slnOj^ 


It  Is  readily  demonstrated  that  the  value  of  Y  In  Figure  11  equals 
the  "effective"  mean  free  path  for  energy  transmission. 

Application  of  these  results  to  dose  measurements  makes  use  of  the 
fact  that,  within  broad  energy  limits,  dose  is  approximately  proportional 
to  energy  flux. 


Appendix  II 

IN- SCATTERING  CORNER  LIP  EFFECT 


Figure  12  llluatrates  the  geometry  considerations  In  the  "In¬ 
scattering  comer  lip  effect  analysis"  (hereinafter  termed  the  "In- 
scatter  effect"). 

The  photons  are  assumed  to  strike  the  comer  In  an  essentially 
parallel  beam,  perpendicular  to  the  lip  height.  This  la  believed  to  be 
a  valid  assumption  for  L/H  ^  3  and  for  a  small  scattering  volioM. 

Let  us  consider  a  photon  striking  at  point  B,  suffering  a  scatter¬ 
ing  collision  at  B  and  Chen  traveling  along  line  BA  to  the  detector. 

If  we  define  a  distance  y  by  letting  AB  a  y,  we  will  concern  ourselves 
for  Che  mament  with  all  other  photons  which  travel  the  same  distance  y 
within  the  comer  lip  material  and  reach  the  detector.  The  other  extreae 
would  be  a  photon  traveling  along  DC  from  Che  source  and  scattering  Just 
°  at  C  Into  Che  detector  such  chat  CD  «  AB  a  y.  The  rays  are  essentially 
parallel  on  entering  or  leaving  the  corner. 

Suppose  we  now  take  any  other  photon  traveling  along  a  path  FG-GE 
which  has  a  scattering  collision  along  the  straight  line  CB  and  Is 
scattered  In  a  direction  parallel  to  the  first  two  photon  exit  direc¬ 
tions.  Will  It  also  have  the  same  distance  of  travdl  In  the  lip 
material?  If  so,  we  can  then  analyse  a  small  volume  element  along  the 
line  CB  and  then  Integrate  over  all  volume  elements  (have  all  path 
lengths)  to  determine  Che  total  single  scattering  contribution  from  the 
comer. 

In  Figure  12  we  are  given  that  AB  ^  CD  ^  y.  The  entering  and  exit 
rays  are  parallel.  The  problem  is  to  prove  that  CD  Is  a  locus  of  all 
scattering  points  which  provide  the  lip  travel  distance  Y  for  all  rays 
having  entrance  and  exit  directions  Indicated. 


MBC  ~;iEGC 

(II-l) 

^iCDB  -  AGBF 

(II- 2) 

Therefore, 

EG  CG 

AB  CB 

(11-3) 

and 

GF  GB 

CD  “  CB 

(11-4) 
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But 


(11-5) 


CB  >  CG  4-  GB 


Therefore , 

EG 

_  (AB)  (CG) 

cr 

(11-6) 

and 

GF 

(CD)  (GB) 

“  CG  +  GK  • 

(11-7) 

But 

AB 

*  CD  -  y 

(II-8) 

Therefore, 

EG  -f  GF 

y  (CG)  .  y  (CB) 

CG  +  GB  CG  +  GB 

(11-9) 

Thus 

EG  4  GF 

-  y 

(11-10) 

It  can  be 

shown  that 

no  point  off  the  line  CB  has 

the  desired 

property.  Therefor*  CB  1*  the  ateted  locus,  which  was  to  be  proved. 

Reconstructing  Figure  12  to  Indicate  the  prlnclpel  quantities,  we 
have  Figure  13,  from  which  the  following  identities  are  obtained; 

OC.-^  (ll-ll) 


y  *ln(B-oi;,) 

00  - 


OB 


y  sina2 
slnR 


(11-13) 


OA 


y  slnCB-Clj) 
sinB 


(11-14) 


CB 


[ 


cotijr  +  cot(&-<f) 


(11-15) 
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or 

CB  =  hx 

(II-16a) 

where 

h  =  ^cotdf  +  cot(8-f)  j 

(II-16b) 

From  AOCB: 

(CB)^  .  (OC)^  +  (OB)^  +  2  (0C> 

(OB)  coa8 

(II.17a) 

.  — Z—  faln^ct  +  atn^O, 
eln^  L  ^ 

n 

(II-17b) 

+  2  slnOj^  alnci2  cosB 

* 

or 

“  alny 

ft 

(II- 18a) 

whe  re 

r  2  2 

J  3  sin  0^  sin  -f  2  alnotj^ 

alnOtj  coscj 

(n-18b) 

Therefore, 

from  Equations  II-16a  and  II>18a, 

(11-19) 

and 

X  h  aln5 
y  =  J 

(11-20) 

or 

y  »  m  X 

(11-21) 

where 

h  8ln& 

“•  j 

(T.I-22) 

since  m  la  a  function  of  the  varloua  duct  angles,  It  can  be  considered 
as  a  constant  In  any  specific  case. 
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with  the  basic  quantities  defined  we  can  now  proceed  to  determine 
the  In-scattering  caused  by  the  corner  lip.  Figure  14  Illustrates  a 
differential  volume  dV,  which  la  equal  to  the  thickness  dx  times  the 
height  of  the  duct  K  times  the  length  CB,  or 


But 


Therefore, 


dV  »  (CB)  H  dx 

(11-23) 

CB  =  h  X 

(II-16a) 

dV  =  H  h  X  dx 

(11-24) 

The  number  of  electrons  par  unit  vt>lume  la  NZ,  where  N  Is  the 
number  of  atoms  per  cc  and  Z  Is  the  number  of  electrons  per  atom.  The 
number  of  scatterers  In  volume  element  dV  •  N  Z  H  h  x  dx,  (11-25) 

If  Iq  la  the  Intensity  at  unit  distance  from  the  source  S,  then 
the  Intensity  at  the  differential  volume  would  be  lo/Ri^,  where  R]_  Is 
the  distance  from  the  source  to  the  scattering  volume,  The  differential 
Intensity  dl  at  the  detector  from  volume  dV  would  then  be; 


I  N  Z  H  h  K(e,E  )  X  e"^®“  dx 
di  .  -2 - S - 

n  2  „  2 

Rj 


(11-26) 


where  e  Is  the  probability  of  a  photon  not  being  absorbed  or 

otherwise  eliminated  from  an  accepted  path 

K(0,E  )  Is  the  probability  of  an  amount  of  energy  being  scat- 
°  tered  through  the  appropriate  angle  Into  the  detector, 
If  initial  photons  have  a  monoenergetlc  value,  of  E 
(Kleln-Nlahlna  Scattering  Probability)  ° 

Rj  Is  the  distance  from  dV  to  the  detector 
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K(0j  .E^) 


(1  +  coa^e) 


2  1  +  7^  (1  -  COS0) 


{ 


(11-27) 


7o^(l  -  .cobS)^ 


(1  +  cos  e)  1  +j'^  (1  -  cose) 


where  Is  the  "dess leal”  radius  of  the  electron 

(r  ^  «  7.9403  x  lO"^^  cm^) 
o 

B  is  the  angle  of  scatter  and  equals  5  -  -  (>2  (II~28) 

y  Is  the  photon  energy  In  terns  of  the  aelf-energy  of  the 
^  electron  and  equals  E^(Hev)/0.511  (11-29) 

For  practical  purposes,  Then  by  Integrating,  we  find; 

^o  o 


H  Z  H  h  K 

°  “  I /*  •* 

^2  o 


dx 


(11-30) 


or  D 


D  N  Z  H  h  K 
o 

„  2  „  i  ,  72 

Rj  Rj  (n») 


(11-31) 


The  volume  which  will  contribute  largely  to  the  In-scattcrlng 
effect  will  be  snail  for  practical  cases.  Figure  15  depicts  the  various 
quantities  required  to  solve  for  D.  fij  and  Oto  sre  the  same  angles  as 
in  the  transmission  comer  effect  (see  Appendix  1),  From  Figure  15, 


2  8 


lnO£^ 


(11-32) 
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(11-33) 


*2  “  i  aina. 


A*  In  prnvloui  easaa,  wc  will  uae  for 


Fron  th«  gacnetry  of  the  duct  it  can  be  shown  thatt 

h  «tna^ 

^2  “  slnB 

Kaklng  this  aubstltutlon.  Equation  3  becoaea: 

Z  N  K  R  sinO^  slnci^ 

D  =  j  ^  j 

4a  *1  *2 


(11-34) 


(11-35) 


With  the  aasuaptlon  that  slnoi  Sr  tanOt,  and  substituting  the  proper  func¬ 
tions  into  Equation  4,  the  final  fora  |or  the  dose  contribution  fron 
in-scattering  in  the  corner  for  S  •  90°  would  be; 


.  *2  «  N 

r  *  **  *  1 

\  2  L,  ) 

.2  4  ^  (1  -  0  )^(1  -  B  )\ 

(11-36) 


Table  IV  contains  values  of  Z  H/4.^  for  various  gama  energies 

and  nuterials.  For  obtaining  values  for  E/r.^,  use  Figure  3.  Assuming 
a  small  scattering  volume,  the  value  of  6  which  should  be  used  for  entry 
into  Figure  3  would  be: 

0^  0  8  -  a,  -  Oj 

(11-28) 

.  -l  ''i 

i  (1  -  Pj) 

(11-37) 

-l  «2 

*^2  ”  2  bj  <1  - 

(11-38) 
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Appendix  III 


FIRST  LEG  ATTENUATION  FACTOR, 

BASED  ON  DOSE  3  FEET  ABOVE  INFINITE  PLANE 


At  the  distance  down  the  first  leg,  the  dose  from  the  rectangular 
entranceway  with  the  Isotropic  source  distributed  equally  across  It  la 
practically  the  same  as  that  from  a  circular  disc  of  the  same  ares,  pro¬ 
vided  the  variation  between  the  values  of  Wj  and  H  Is  not  too  extreme. 

A  disc  of  this  area  would  have  a  radius  equal  to  H/a.  From 

Section  4.12  of  Reference  6,  Equation  4.12,1a  Indicates  that  the  dose 
at  distance  L  Is  approximately: 


D 


1 


,  H 
2ir  L^^ 


(III-l) 


where  J  la  a  factor,  depending  on  source  strength  and  energy,  which 
provides  numerical  answers  In  proper  units  of  dose  or  dose  rate. 

A  rather  accurate  formula  for  dose  readings  In  air,  3  feet  above  a 
smooth,  Infinite  plane,  homogeneously  contaminated  by  a  distributed 
radiation  source.  Is  given  (Ref.  9)  as  follows; 


I  r 

r  - 


-3Bi^^alr  .  . 
e  +  A„ 


-3B2Walrl 


(I1I-2) 


where  Ugir  Is  the  total  absorption  coefficient  in  reciprocal  feet  for 
the  radiation  in  air;  A2.  ,  and  B2  are  empirical  coefficients  to 

account  for  the  scattering  contribution  (Ref.  10);  and  Is  the  ex¬ 
ponential  integral  (Ref.  11). 

Using  an  energy  value  of  1  Mev  for  the  photons,  we  can  readily 
calculate  that 

D  a  5.36  J 
P 


(III-3) 


The  factor  5.36  la  dependent  on  photon  energy,  but  to  a  rather  slowly 
varying  extent.  It  will  therefore  be  sufficiently  accurate  to  consider 
it  as  representative  for  any  practical  situation  of  this  sort. 

Then  the  first  leg  attenuation  factor,  baaed  on  the  dose  at  a 
reference  point  3  feet  above  a  smooth,  Inflnlt#,  contaminated  plane,  is; 
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(III-4a) 
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